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y-^ , Higher twist effects in the deeply inelastic scattering are studied. We start with a short review of 

the theoretical results on higher twists in QCD. Within the saturation model we perform a twist 
analysis of the nucleon structure functions Ft and Fl at small value of the Bjorken variable x. 
The parameters of the model are fitted to the HERA F2 data, and we derive a prediction for the 
longitudinal structure function Fl- We conclude that for Fl the higher twist corrections are sizable 
whereas for F2 = Ft + Fl there is a nearly complete cancellation of twist-4 corrections in Ft and 
'""^ ' Fl- We discuss a few consequences for future LHC measurements. 
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I. INTRODUCTION 

A deeper understanding of the transition region at low Q 2 and small x in deep inelastic electron proton scattering 
■ has been one of the central tasks of HERA physics. Approaching the transition region from the perturbative side 
£f~) [ one expects to see the onset of corrections to the successful DGLAP description — in particular those which belong 
0^ ■ to higher twist operators in QCD. The twist expansion defines a systematic approach and, therefore, provides an 
attractive framework of investigating the region of validity of the leading twist DGLAP evolution equations. 

The essentials of the theory of higher twist operators and their Q 2 evolution have been laid down several years ago, 
U S HI- First, a choice has to be made of a complete operator basis @, and for their evolution [l[ one needs to 
compute evolution kernels which, for partonic operators in leading order, reduce to 2 — ► 2 scattering kernels. The 
problem of mixing between different operators has also been addressed first in [l[ . In the small- x region at HERA 
one expects the gluonic operators to be the most dominant ones; so far, a theoretical study of the evolution of twist-4 
gluon operators is available only in the double logarithmic approximation 0, S S 0| ■ An extensive recent theoretical 
study of QCD evolution of the higher twist operators can be found in Q . 

Numerical studies of the size of potential higher twist corrections 0, licill indicate that twist-4 corrections to Fi 
are small down to Q 2 ~ 1 GeV 2 , x ~ 1CP 4 . A first theoretical analysis [llj applied to HERA data, however, has 
shown that the situation is more subtle, and from the smallness of twist-4 corrections to F2 one cannot conclude that 
contributions of twist-4 operators are small. The simplest QCD diagrams contributing to the twist four gluon operator 
are shown in Fig.[TJ a quark loop couples, via the exchange of four i-channel gluons, to the proton target. Calculating 
the contribution of these diagrams to the 7*p cross section at small x, and isolating the twist-4 contribution one 
arrives at the conclusion that the contributions to longitudinal polarized photon has the opposite sign compared to 
the transversely polarized photon. This implies the possibility that, in F2 = Fl + Ft which sums over transverse and 
longitudinal photons, there is a (partial) cancellation of twist four corrections, whereas the twist-4 corrections to Fl 
or Ft are larger than the corresponding corrections to F2. 

In order to decide whether the HERA data support this possibility, in [ll[ the saturation model of Golec-Biernat 
and Wiisthoff [l2|, HH which successfully describes the HERA data has been used to obtain a quantitative estimate 
of twist-4 and even higher twist corrections. This mathematically fairly simple model has four parameters which 
are fixed by adjusting the model to describe well the HERA F2 data. The model leads to the total cross sections 
that exhibit the geometric scaling (T3 |. The cross sections obtained in the model may be expanded in powers of 

^ (j) A - It was natural to identify the first two terms of this expansion as 'leading twist-2' and 'twist-4 correction', 
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respectively. Despite its simplicity this model is close enough to the lowest order QCD calculations and supports 
the sign structure of the twist-4 correction to the transverse and longitudinal cross sections mentioned before. On a 
quantitative level, the twist-4 corrections to Fl and Ft were found to be sizable, whereas in F2 they almost cancel. 
The overall smallness to the 1 /Q 2 corrections to F2 is consistent with the estimate of d, |T(| : the analysis in [ll[ 
therefore provides a natural explanation of the suppression of twist-4 corrections to F2 , without demanding that the 
higher twist contribution is small for Fl or Ft- 

This original version of the GBW model did not include any QCD evolution. Therefore, the connection of this 
model with evolution of twist four operators in QCD is not possible, and, in particular, the identification of the first 
and the second term in the 1/Q 2 expansion as the 'leading' and the 'next-to-leading' twist seems somewhat crude. 
A more recent version [15( of the GBW model includes QCD evolution and its description of HERA data is slightly 
better than that of the original model. It is therefore natural to investigate to what extent this improved model 
exhibits the structure expected for higher twist operators, and then to perform a numerical analysis similar to [Til ] . 
This is the goal the present paper. 

Our numerical analysis shows an interesting pattern of the higher twist effects in the structure functions. The 
corrections are sizable in Fl, for the kinematic range relevant for HERA data at low Q 2 , 1.8 GeV 2 < 10 GeV 2 
and small x, the twist-4 corrections are found to reduce the leading twist result by about 20-50%. We compare 
the obtained predictions for Fl to recent HERA data both for the complete saturation model and its leading 
twist component. For F2 the higher twist effects are found to be surprisingly small, at a few percent level down to 
Q 2 = 1 GeV 2 . 

The paper is organized as follows. The first, longer part in Sees. II- VI is devoted to discussion of theoretical 
issues, and the second part, in Sees. VII-IX , to the phenomenological applications. In section II we review simple 
QCD calculations of higher twist corrections in momentum space, restricting ourselves to the double logarithmic 
approximation, and in section III we reformulate these results in the QCD dipole picture. Next we turn to the 
saturation model: after a brief review of the simple GBW model in section IV, we perform a theoretical twist analysis 
of the QCD improved dipole model in section V-VI. The second part (section VII-IX) contains our numerical analysis 
and discussions of the results and possible consequences for physics at the LHC. Conclusions are given in section X. 



II. THE DOUBLE LEADING-LOG APPROXIMATION IN QCD 

In this section we give a brief overview of higher twist corrections in small- a; QCD. We consider the scattering 
of a virtual photon with transverse or longitudinal polarization on a quark, and we restrict ourselves to the leading 
behavior at large Q 2 and small x (double logarithmic approximation, DLA). In this limit, we can either start from the 
leading-log Q 2 limit and then take the small- 2 limit or, alternatively, start from the small- x limit and then investigate 
the large-Q 2 approximation. We begin with the latter one, i.e. we restrict ourselves to those diagrams which have the 
maximal number of logarithms in X/x. 



A. Twist 4 corrections 

In the small- 2; limit, the scattering of a virtual photon on a quark is described by the exchange of 2 or more gluons 
between a closed quark loop and the target quark. For two ^-channel gluons, the coupling to the virtual photon 
is described by the well-known photon impact factor, -D(2,o);T.l(Q; &j Q ~ k), where the subscripts T, L refer to the 
polarizations of the virtual photon, and k and q — k are the transverse momenta of the two i-channel gluons. In the 
deep inelastic limit: q 2 = 0, k 2 <c Q 2 we have [17l[iq : 

D im;T (Q;k,-k)= a ( f^iogg + K^ )2 +o{{*)*j\ 

where Aq = \^2e 2 a s ^. The leading power, (fc 2 /Q 2 ), belongs to leading twist, the terms proportional to (fc 2 /Q 2 ) 2 
to twist four etc. These results, obtained directly from Feynman diagrams in the momentum space, lead to estimates 
of twist contributions to Ft and Fl consistent with estimates of the saturation model, discussed later, in Sec. IIIII 
and Sec. IIVI obtained in the Mellin representation. Important features of -D(2,o) -t,l are the logarithms and signs 

of the (fc 2 /Q 2 ) 2 corrections: whereas for transverse polarization there is no logarithmic enhancement and the power 
correction is positive, for longitudinal polarization we have a logarithmic enhancement, and compared to the transverse 
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FIG. 1: Two diagrams with four f-channel gluons. 



case is has the opposite sign. These simple observations open the possibility that, in Fl, the higher twist corrections 
are large and that, in the structure function F 2 which sums over transverse and longitudinal polarizations of the 
photon: 



F) — Ft 



Ft 



(2) 



the total twist-four corrections may be small due to cancellations. Within the model to be discussed in this paper we 
will find that, for the HERA data, this is indeed the case. 

The simplest diagram for corrections due to the exchange four gluons (order 0(g 8 )) is illustrated in Fig.[TJ left. In 
addition to this exchange diagrams, there are diagrams involving the triple gluon vertex, like e.g. the diagram shown 
in Fig. [H right. 

An efficient method of calculating the sum of all these contributions in the high energy or small- a; limit is the 
effective action, defined in [III] and further studied in [2(|. As a result we have, for all diagrams up to order g 8 , two 
classes of contributions: 

(a) the BFKL ladders [21, 22, 23, 24], expanded up to order g 8 . This includes NLO corrections to the gluon trajectory 
as well as NLO corrections to the BFKL kernel which, in the leading log approximation, will not be considered; (b) 
the exchange of four i-channel gluons, where the eikonal couplings to the quark loop at the top and to the quark line 
at the bottom arc fully symmetrized, both in color and in their momentum structure. The coupling to the quark 
loop, -D(4^). T l^i, k%, ^3; ^4), can be expressed in terms of D 2 : 



%,0);T,LX 2a3a4 (fcl,fc2,fc 3 ,fe4 



^(2,0);T,i(*!l + fc 2, + &3 + ^4) 



. , _ 2 101020304 
"4 1 — g u sym 

(-D(2,0);T,i(fel + k 2 + k 3 ,k 4 ) + £ ) (2,0):T,L(fc2 + &3 + ki,kl) 

+ £>(2,0);T,L(fc3 + fc 4 + fcl, fca) + £>(2,0);T,L (&4 + &1 + &2, fes) 

^(2,0);T,i(fel + fe 3, &2 + ^4) - D(2,0);T,L (&1 + &4, &2 + £3)) , 



(3) 



where <ij and ki , i = 1 , . 
form: 



,4 are the color indexes and transverse momenta of the gluons. The color factor has the 



101020304 

* sym — 



tr (t ai t a H a n ai ) + tr (t a H a H a2 t ai ) 



(4) 



where the subscript 'sym' indicates that the color labels are completely symmetrized. 

Next let us consider higher order corrections to (a) which, in the leading logarithmic approximation, sum up to 
the BFKL ladders. Using a Mellin representation for the impact factors D/ 2 ,o);T,L(Qlk,—k) and a double Mellin 
representation for the BFKL Green's function we write the scattering amplitude in the form: 



W^7QS> = /^/£(iH|) W.M 



1 



(5) 



where the integration contours run along the imaginary axis, a = ^ 
form 



U3 - ax{v,0) 
the BFKL characteristic function has the 



x(v,0) = 2il>Q)-ili(l + v)-iK-v), 



(6) 



and Qq denotes the momentum scale at the target end of the BFKL ladder. The scattering amplitude D 2 can be 
expanded in powers of Q 2 /Qq'. the terms in this expansion are due to the poles of xiy) a ^ negative integer values: 
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FIG. 2: (a) Pairwise interactions between four reggeized gluons; (b) two noninteracting color singlet ladders. 

v = — 1, —2, and the residues of the poles lead to anomalous dimensions. In particular, the pole near v = — 1 leads 
to the leading twist behavior 

[Q /Qo) > 7M = (7) 
the pole near ^ = — 2 to the twist four correction 



TTUO 



[Q /Qo) , tM = (8) 

The exponents are the anomalous dimension of the higher twist operators contained in the BFKL amplitude [25[ . 
In the notation of they belong to the class of non-quasipartonic operators. When coupled to the impact factors 
D(2,o);T,l{Q] — k), we see from the expansion in (TTJ) that, for the leading-twist terms, the transverse polarization has 
a logarithmic enhancement compared to the longitudinal polarization. For the twist four correction, the situation is 
reversed. As to the sign structure, the twist four corrections to transverse and longitudinal polarizations have opposite 
signs. 

Next we turn to higher corrections to class (b) . Following the analysis of [l[ we consider those corrections which are 
obtained by inserting all possible pairwise interactions between the reggeized i-channel gluons (Fig. [5^). The large-Q 2 
behavior comes from the region where the transverse momenta of the exchanged gluons are ordered and much smaller 
than the virtuality of the external photon. The leading (in 1/Q 2 ) behavior of the four gluon state has been discussed 
in 0,0,0], and we briefly summarize. In the complex z^-plane, the leading singularity is a pole at 

_ + 
iruj 

where S is a correction of the order (5 = 0.778/N*, and for N c = 3, S = 0.0096). In the large- iV c limit the four 

gluon state reduces to two noninteracting color singlet gluon ladders, leading to a cut in the z^-plane with the branch 
point located at 

lent = 10 

TTUJ 

The pole in © at finite iV c can be viewed as 'bound state' formed by the two color singlet ladders, whereas the cut 
(HHJ) represents the 'threshold' of two free ladders. The large-Q 2 behavior of this four gluon state is described the 
evolution equations of the twist four gluon operator in the small- a; limit, as discussed in In the large- iV c limit, 
the evolution equations reduce to two independent DGLAP ladders. 

In order to apply this discussion to the diagrams of class (b) we notice that the twist four contribution of the 
fermion loop, Du \ T L is easily obtained from ([3]) and Jl]). For example, the twist four correction of the transverse 
polarization is found to be proportional to 



2 2 ki ■ k 3 k 2 ■ &4 + fci • k 2 k 3 ■ fc 4 + ki ■ fc 4 fc 2 • k 3 
49 Ao 5 W^ 2 ' ( } 



and analogous results are found for the longitudinal polarization. In this way, the four gluon state corrections to 
transverse and longitudinal polarizations follow the same pattern as the twist four piece inside the BFKL ladder: 
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compared to transverse polarization, the longitudinal polarization has a logarithmic enhancement and comes with the 
opposite sign. 

Compared to the BFKL-singularity in 1(5)), both in © and in (fTUf the coefficients of the pole at co = are larger 
by a factor 4: at small x, this twist-4 correction will therefore dominate. This suggests to consider, within a twist 
expansion in the small- a; region, as a first set of higher twist corrections these four gluon states, disregarding the 
higher twist contributions of non-quasipartonic operators. It is not difficult to generalize this selection to six, eight 
etc gluon states. If, in addition, one invokes the large- N c expansion where the 2n-gluon state is approximated by n 
noninteracting color singlet ladders, one arrives at the eikonal picture of multi-Pomeron exchange, which underlies 
the saturation model to be discussed further below. We shall see that this model embodies many of the features of 
these n-ladder exchanges, in particular the correct Q 2 -evolution. 



So far we have discussed a selected subclass of QCD diagrams giving rise to twist-4 corrections to the proton 
structure functions. It should, however, be kept in mind that this selection of higher twist corrections is not in agree- 
ment with what one obtains from summing all leading-log 1/x contributions from the BFKL Pomeron fan diagrams. 
This summation may be performed using the Balitsky-Kovchegov (BK) equation [2S [27|, 28]. To illustrate this, we 
summarize the results of the complete small- 2; analysis which can be found in [il. l29l |30|. The sum of all diagrams 
contributing to the leading logarithmic l/x approximation can be organized in two classes. 

(i) BFKL ladders consisting of reggeized gluons. At the lower end, reggeized gluons can split into two or three ele- 
mentary gluons. In the former case, the color structure of the splitting is described by a structure constant f cai a 2 i m 
the latter case by the product of two structure constants, e.g. f C a 1 dfda 2 a 3 - 

(ii) BFKL-like ladders where, instead of the reggeized gluons which belong to the adjoint color representation, we 
have Reggeons in the symmetric octet and singlet color representations. In both cases, the trajectory functions are 
the same as for the reggeized gluon. The corresponding color tensors are listed in [2i|. The sum of these diagrams 
is symmetric under the exchange of momenta and color indexes. 

(iii) Diagrams with a four gluon t-channel state. This state is symmetric under the exchange of t-channel gluons 
(momenta and color indexes). There is no direct coupling of this state to the quark loop at the top. Instead, through 
the 2^4 reggeized gluon vertex it couples to a BFKL ladder which then connects with the quark loop (Fig. [3]) . Class 
(i) and (ii) represent the all-order generalizations of (a) and (b), respectively, whereas (iii) starts at the order g 10 . In 
[2g |. (i) and (ii) are denoted by Df, class (iii) by D\. 

As we have already stated, in (iii) the four gluon state that we have discussed before does not couple directly 
to the quark loop: the coupling goes through a BFKL ladder and a 2 — » 4 transition vertex. Making use of the 
large-Q 2 results discussed before, we interpret this as a mixing between the non-quasipartonic twist four piece inside 
the BFKL ladder and the twist four gluon operator. A detailed analysis [3(j, however, shows that, at the leading 
logarithmic logQ 2 , approximation, this transition kernel between the two twist-4 operators, in the large- N c limit, 
vanishes. This also holds for the transition of the twist-6 piece inside the BFKL ladder to the twist-6 piece in the four 
gluon state. Generalizing this to more than four t-channel gluons, one arrives at the conclusion that, in the double 
leading logarithmic approximation, the contribution of higher twists given by the BK fan diagrams vanishes, and only 
propagation of two interacting t-channel gluons contributes to the amplitudes. Note however that this result is only 
valid in the large- N c limit. 

Therefore, from a theoretical point of view, the twist expansion derived from the selection of diagrams with 2n 
t-channel gluons in the large-7V c limit should be viewed as being different from descriptions based upon the BK 
equation. 



It is important to emphasize that, in the small- a; limit, for the class of QCD diagrams which we have discussed, the 
scattering amplitude for the elastic scattering of a virtual photon on a quark can be cast into the dipole form [3l| : 



B. Higher twists in the Balitsky-Kovchegov equation 



III. COLOR DIPOLE PICTURE 




(12) 
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FIG. 3: A fan diagram. 



where T and L denote the virtual photon polarization: transverse and longitudinal, respectively. The light-cone 
photon wave function, ^ T L , is modeled by the lowest order 7*5 — > qq scattering amplitudes which give 

\H(z,r,Cf)\* = 2N 2f" { [* 2 + (1 - *?] ^ Kt{er) + mj *g(er)} (13) 
|*{for,Q 2 )| 2 = ^OV(l-f^H (14) 

where -Ko,i are the Bessel-McDonald functions, e 2 = z(l — z)Q 2 + raj and r = |r|. The measured structure functions 
are related to crJ£(x,Q 2 ) by the standard formula 

FT - L = (15) 

In (|12p , all details describing the interaction of the quark-antiquark pair with the target quark are contained in the 
dipole cross section, <j(x, r). In particular, the exchange of two non-interacting color singlet gluon ladders provides a 
contribution proportional to the product of two gluon structure functions, (xg(x, C/r 2 )) 2 . 

Important characteristics of the twist expansion follow from the structure of the photon wave functions and do not 
depend upon the details of <r(x,r). This is most easily seen by taking the Mellin transform of (|12[) . In general, the 
Mcllin transform of a function /(r 2 ) is defined as 

/>oc 

f( S )^M r2 [f(r 2 )]( S ) = / dr 2 {r 2 y- 1 f{r 2 ) (16) 

Jo 

while the inverse relation reads 

f(r 2 ) = j c ±L{r 2 )-°f{s) (17) 

where the integration contour C lays in the fundamental strip of the Mellin transform to be discussed below. 
Let us write Eq. (fT2|) in the following form 

,-r" ''. ' 

where 



all(x,Q 2 ) = I —H T , L (r,Q')a(x,r) (18) 



H T , L (r,Q 2 ) = irr 2 V f dz\¥ T L (z,r,Q 2 
t Jo 



(19) 



/ 

Substituting the inverse Mellin transform of the dipole cross section, 



"for) = J ±,(r 2 Q 2 )- s a(x,s) (20) 
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we find the Mellin representation of the j*p cross sections, given by the Parseval formula 

4l{x,Q 2 ) = f ^ta(x,s)H T , L (-s,Q 2 /Ql) (21) 

J 2lTl 

c. 

where Ht,l{—s,Q 2 ) is the Mellin transform of HT,L(r,Q 2 ). The integration contour C s in the complex s-plane is 
placed in the fundamental strip in which the integrals defining a(x, s) and Ht,l{—s, Q 2 ) are convergent. The strip is 
determined from the following leading behaviour of both functions at small and large values of r (up to logarithms of 
r): 

tt / ^>2s f const for r — > 0, /r<r , s 

HT,L(r,Q 2 ) = { 1/r 2n for ^ (22) 

with n = 1 for transverse and n — 2 for longitudinal polarization. For the dipole cross section we take, as an example, 

I \ ! r2 for r 0, . . 

a(x,r) — < , t (23) 
v ' I const for r — > oo v ' 

In this case the fundamental strip of a(x,s) is defined by the condition — 1 < Res < while the fundamental strip 
of Ht,l{ s iQ 2 ) is given by < Res < n. Taking into account the minus sign in Ht,l(~s,Q 2 ), we find that the 
integration contour C s in Eq. (f2"Tj) lays in the strip: 

-KReKfl. (24) 

It can be chosen parallel to the imaginary axis, for example, s = — 1/2 + iv with real v. 

For the Mellin transform of HT,L(r,Q 2 ) we restrict ourselves to massless quarks, m,f = 0. In this case, Ht.l are 
functions of only one combined variable, f = rQ: 

H T (f) = A [ dz [z 2 + (1 - z) 2 ] z(l - z) f 2 K 2 (y/z(l -z)r) (25) 



H L (f) = 4A / dzz 2 (l-z) 2 r 2 K 2 (^z(l-z)r). (26) 



where we introduced A$ = N c a em (e 2 } / (2ir) and (e 2 ) = X)/ e /- ^ n Appendix [Al we found the Mellin transforms of 
these functions in the form 



#t,l( S ,Q 2 ) = ( ? ) H t,l{s) (27) 
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with Ht,l(s) given by Eqs. (TA"9j) and ]HD| : 



and 



- Aott r(2 + s)r(i + s)r(s)r(i-s)r(3-s) 

Tl ' J 8 r(3/2 + s)r(2-s)r(5/2-s) ' 1 ' 



h m A ° n (r(i + ^)) 3 r(2-s) 

Hl{s) - — r(3/2 + s)r( 5 /2-s)- (29) 



Both functions have simple or multiple poles for negative and positive real values of s. 
Substituting these results into Eq. (l2Tj) we obtain 



aU(x,Q 2 ) = I £ (%) ~a{x,s) H T , L (-s) (30) 



2iri \Q 2 



with the contour C s in the fundamental strip (j2~4"|) . The twist expansion is obtained by closing the s-contour to the 
left. The functions Ht,l{s) have single poles to the right of C at positive integers, except for the regular points 
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at s — 2 for transverse and s — 1 for longitudinal polarizations. Thus, both functions have the following Laurent 
expansion around each singular point s = n: 

(n) 

H T<L (s) = + bfrl + O(s-n) (31) 



with 



4 2 ) = 0, = --A (32) 



and 



0, b^> = \Ao 



«i 2) = -g^o. (33) 



Assuming, for simplicity, that to the left of C s , a(x,s) has only poles at s = —1,2, ... (the more realistic case where 
a has cuts in the complex s-plane will be discussed further below), we close the contour to the left and arrive at the 
twist expansion: 



oo 



, t m 2 ) = £4:r 2n) M 2 )- (34) 

n=l 



where a^ T=2n ^ ~ l/Q 2n (modulo powers of logQ 2 ). With {25} and we obtain: 



r (r=2„) _ f ds ( Ql\ s ,^ ^ J -Ot"l ^_ fc ( tt ) 
27ri VQ 2 



c„ 

where the dots stand for terms regular at s = —n. In particular, for the twist-4 corrections we re-discover the previous 
result from ([T]): 

(2) 

(i) due to the vanishing of a T , the longitudinal structure function is enhanced, 

(ii) the leading terms in Ft and Fl come with opposite signs. 

For completeness, we also consider the complex half s-plane to the right of the contour C. It is well known that 
the Bessel-McDonald functions K v (x) have a convergent expansion around x = 0, whereas for large arguments the 
expansion in powers of 1/x is asymptotic. Therefore, writing the functions Ht,l{^) in the form 

_ A n r d s (? 2 y s r(2 + a )r(i + a )r(s)r(i- s )r(3- 8 ) 

T[ ' 8 J2iri\iJ r(3/2 + s)r(2-s)r(5/2-s) 1 ' 

c 

Aqtt r d s (py s (r(i + .)) 3 r(2- 3 ) 

Hi(r) = — i2^UJ r(3/2 + .s ) r(5/2- S ) ' (37) 
c 

we conclude that the expansion in powers of f — which is obtained by closing the contour to the left — is convergent. 
In contrast, the expansion in powers of 1/f — which corresponds to closing the contour to the right and computing 
residues of the poles at positive integers — leads to a divergent result which form an asymptotic series for Ht,l{?) 
when f 2 — * oo: 



h {1) h {2) h (3) 

H t ,l{t) ~ 4^ + ^ + ^ +.-., (38) 



where the coefficients ~ a TL are ec l ua l to: 



4 X) = | A) , = » h T ] = W A o (39) 

fc«=0, 4 2) = W^o, ^ 3) = ^o. (40) 
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This asymptotic expansion justifies the large-r behaviour of Ht.l (r, Q 2 ) used in the determination of the fundamental 
strip (|M|) . Moreover, returning to (f2"Tj) we conclude that, because of the negative sign of the argument of Ht,l, the 
twist expansion is an asymptotic expansion. 

In conclusion, the opposite sign structure as well as the relative enhancement of the twist-4 corrections to Fl are 
general features of the small- a: limit in QCD, and they provide the possibility that the total twist-4 correction to F 2 
may become small. In the following we choose, for a quantitative estimate, a particular model, the QCD improved 
dipole model. 



IV. THE MODEL 



We aim for the construction of the twist expansion of the proton structure functions Ft and Fl at small values of 
the Bjorken variable x. The starting point for our following analysis is the GBW saturation model [l2j and its QCD 
improved version which incorporates the leading logarithmic DGLAP evolution [l5j |. 

The standard formula defining the total cross section for the scattering of a virtual photon j T L (Q 2 ) on a proton p 
at small value of the Bjorken variable x has already been written down in (fT2|) . The function <j(x, r) in Eq. (fl2|) is the 
color dipole cross section, describing the interaction of the qq pair with the proton. In the original GBW formulation 
[l2| it depends on the dipole size r and the Bjorken variable x, and takes the following form 

a(x,r) = a Q {l - exp (-r 2 Q s 2 at (x)/4) } (41) 

where Q 2 at is a saturation scale which depends on x. After incorporating the DGLAP evolution for small dipole sizes 
the dipole cross section is modeled in [15| as 

a(x,r) = cr {l-exp(-O(x,r 2 ))} (42) 

where the opacity 

n( x y) = ^ 2 ^(M 2 )g(*,M 2 ) (43) 

3(T 

and g(x, /1 2 ) = xG(x, fi 2 ) is the gluon distribution (multiplied by x) which obeys the DGLAP evolution equation (|B 1 1) 
from Appendix [BJ The evolution scale /i 2 was originally assumed to depend on the dipole size in the following way: 

fi 2 = C/r 2 +p 2 . (44) 

Both models of the dipole cross section are eikonal and follow the Glauber-Mueller formulae. For the remainder of 
this section, we restrict ourselves to the original GBW model. 

Following our discussion of the previous section, we need the Mellin transform of the dipole cross section. In the 
case of the GBW parameterization (|4Tj) . we find 

a(x, s) 



a I dr 2 (r 2 )- 1 {1- exp(-r 2 g 2 at ( 2 ;)/4)} 

[ C AA (45) 



4 

Substituting this result, together with relation l|27[). into Eq. (l2Tj) , we obtain 

af}(x,Q 2 ) = a [ ^ (%) ' {-T(s)} H T , L {-s) (46) 



2ni V Q 2 



with the contour C s in the fundamental strip (j2"4")) . We see that the poles to the left of C s at negative integers lead 
to the twist expansion: 



rl}(x,Q 2 ) = J24= 2n \x,Q 2 ). (47) 

71=1 



10 



where cr( T=2n ) ~ l/Q 2n . Singularities come from the single poles of the Euler gamma function T(s) and from the 
poles in Ht l(— s). In particular, encircling the pole at s = —n by a small counter-clockwise oriented contour C n , 
and expanding both functions around this point, we obtain 

(r=2n) f ds f Q s 2 at y S f 7i" } M \f , (») 1 r4R s 

^ - ao j2ri{-&) \~ + 7 ° \7T^ + Vl + -7 (48) 

where the dots denote terms regular at s = — n. The result of the integration is indeed proportional to 1/Q 2n with 
the logarithmic enhancement coming from the double poles 

4T l) = || {-7i n) 41 log(Q 2 /QL) + (7i n) &S - 7^41) } • (49) 
In particular, we find [ll[ _ for twist-2: 



2 



(T 



(r=2) _ "em^O / 2 \ Wsat 



)2 



( e 2) ^ {l 0g(Q 2 /Q 2 at) + ^ + 1/6} (5Q) 



and for twist-4: 



(t=2) _ ^em^O / 2\ Vsat 



(r=4) _ 3 a em ffo ; 2\ Qsat /r \ 



i T=4) = -^^l( e 2 )%{log(Q 2 /QL)+ 7^ + 1/15}. (53) 



Notice the negative sign of ^ and the lack of logarithm in <j£ ^ and <r^T 4 - > due to the singularity structure (|3Tj) 
with = and = 0. 



V. SINGULARITY STRUCTURE OF THE DGLAP IMPROVED MODEL 



In the DGLAP improved saturation model, the r-dependence of the dipole cross section, given by Eq. R2[) . is rather 
involved and its exact Mellin transform is not known. However, it is still possible to extract the information about the 
Mcllin transform necessary to carry out the twist analysis. For this purpose it is convenient to use a slightly modified 
definition of the scale /i 2 in Eq. 

2 _ f C/r 2 for r 2 < C/ M 2 , 
' ~ \ for r 2 > C/fi 2 . (&4j 

Such a modification preserves all the desired features of the original model and allows to separate the reintegration 
range of the Mellin transform of the dipole cross section B{x, s) into two regions: the perturbative one, defined by the 
condition r 2 < C//Zq, in which the gluon density and strong coupling constant are given by one- loop expressions with 
the scale fi 2 = C/r 2 , and the soft region, defined by the condition r 2 > C//j,q, where the scale is frozen at y 2 = fi^: 
Thus 

a(x,s) = 5- pert (x,s)+a so ft(x,s). (55) 

In the soft region the dipole cross section takes the form of the GBW saturation model ijlTj) with the saturation 
scale 

QL(*) = 4?r2as(/ f g(x ' M§) - (56) 
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FIG. 4: The singularity structure of a(x, s)Ht(-s) in the complex s-plane to the left of the Mellin integration contour C s in 
the fundamental strip shown as the gray band. The zigzag lines indicate the cuts coming from powers of the Mellin transform 
of a s (/J, 2 )g{x, /x 2 ) while the full circles are the poles of Ht(—s) that coincide with branch points of a(x, s). The crossed circle 
at s — —2 is the branch point that is not accompanied by the pole of Ht( — s). 



The contribution from this region to the Mellin transform is given by 

/>oo 

a^x, s) = / dr 2 (r 2 )*- 1 {l - exp (-r 2 Q s 2 at (*)/4)} 
JC/nl 



= - (T0 { 9 t) |ir +r(s ' a) j (57) 

where a = CQ 2 at /(4/ig) and T(s,a) is the incomplete gamma function which has no singularities in the complex 
s-plane. The soft part has only a single pole at s — which does not contribute to the twist expansion. 

The contribution from the perturbative region may be computed term by term from the expansion 



a pcrt (x,r) = ^cr^ t (>,r) 



71=1 



°° (-^Y l + 1 

'£a { —L—n'Z elt (xy) (58) 



n=l 



and the perturbative part of the opacity reads 

tt P crt(x, r 2 ) = ^- r 2 a s {C/r 2 )g(x, C/r 2 ) 9 (C/r 2 - fi 2 ) . (59) 

The Mellin transform <j peT t(x, s) exists due to the theta distribution and is given by the sum of the Mellin transforms 
of the subsequent terms in the expansion 



M (-1)" +1 ~ 
&pelt( X > s ) = — ^pcrt(a;, s) (60) 

Each term contributes a cut singularity in the s-plane extending to the left from the branch point at negative integers, 
see Fig. The positions of the branch points are determined by the corresponding power of r 2 since the powers of 
a s (/j, 2 )g(x, [i 2 ) do not introduce any additional shift. 
For example, we compute first the Mellin transform 

M r , [a s (C/r 2 )g( Xl C/r 2 )Q(C/r 2 - fi 2 )} (s) = (JL\ a^(x, s) (61) 

where a^g denotes the Mellin transform with respect to the scale /i 2 , defined as 

arf(x,s) = [ ^ (j-) a s (ii 2 )g(x^i 2 ) . (62) 
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a s (n 2 )g(x,fi 2 ) = I ( Yn ) « s ff(a;,s). ((>:-!) 



In such a case the inverse relation reads 

^ 2 ).9<>,M 2 ) = / 1 v> 

where the integration contour lays to the right of the right-most singularity. Using the property of the Mellin transform: 
M t [t n f{t)]{s) = M[/(i)](s+n), wefind 

M r 2[n pcrt (x,r 2 )]{s) = — f — J 5^(a;,s + l) (64) 
and after substituting solution (|B16p . we obtain 

<W(*,.) = f - (^) S+1 / ^ ( S + l)-^H. (65) 

3(7 V A / J 2m P gg{u) 

The logarithmic cut singularity along the negative real axis with the branch point at s = — 1 is obvious from this 

} 2 ( 

'pcrtl 

s+2 

{a s g) 2 (x,s + 2) (66) 



solution. The Mellin transform of fl 2 t (x, r 2 ) is given by 



2 \ 2 
— ) 


(- 


3<7 J 






(~ 


3<7 / 


\A 2 



7T 2 \ 2 fCV +2 f ds' . , „ 

3^J [l?) / S M(v)^(M + 2- S ). 



(67) 



where (a s g) 2 is the Mellin transform (|62p of the product [a s (/i 2 )g(a;, // 2 )] 2 , and the Mellin convolution theorem was 
used in the last equality. It can be shown explicitly that expression (1661) has a cut singularity along the real axis for 
— oo < s < — 2 with the branch point at s = — 2. 

In general, we have 

IT 2 \" f C ^ S+n 



MAW; cit (x,r 2 )](s) = I — 1 1^1 (a s g)"(x,s + n) (68) 

with the logarithmic cut along the negative real axis starting at the branch point at s = —n, see Fig. 01 In summary, the 
singularity structure of the Mellin transform (|55p , relevant for the twist expansion, is determined by the perturbative 
part only. 



VI. TWIST DECOMPOSITION IN THE DGLAP IMPROVED MODEL 



At each twist the saturation model incorporates a few distinct contributions that have a clear interpretation 
within perturbative QCD. The contributions may be classified using the singularity structure of the product 
a(x, s)Ht,l(—s, Q 2 ) in the Mellin plane. 



A. Twist-2 contributions 



Starting from the twist-2 analysis, we close the contour C s of the Mellin integration in Eq. (|2"Tj) with two large 
quarter-circles Q\ and Qi and a contour T> s enveloping the complex cut of B{x, s) with the branch point at s = — 1, 

see Fig. [5] Then, we decompose a(x, s) into a part which singular at s = —1, given by cf^ rt (:r, s) = cro^ p0 rt(^, s), 
and a part which is regular at this point, areg ^{x, s). The latter part consists both the soft contribution (|57p and 
the contributions from multiple exchanges with cuts starting from s = —2. Thus, using expansion (|3ip for Ht,l(~s) 
with n = 1, we obtain the twist-2 part in the form 

^ = / ^(^)\^ pcrt (x,s) + ^\ (69) 
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FIG. 5: The discontinuity structure of <7p^ rt (x, s)Ht{—s) and the integration contours in the complex s-plane with the pieces: 
Si, 0.2 and D s . The meaning of all symbols is as in Fig. [4] 

where the integration contour is reversed with respect to the contour D s shown in Fig. [5) and the Laurent expansion 
coefficients 



- 2 A 



a« = 0, 



- 2 A 



(70) 



with A =N c a em {e 2 )/(2Tr). 

Let us compute the twist-2 contribution for transverse photons coming from the most singular part of the Mcllin 
integrand: 



Aa4 T=2) 



2ni V 4 



(71) 



-D, 



The analogous longitudinal contribution vanishes since — 0. Using relation (|64[) . we find 



4m 3 f ds_ (CQ 2 y +1 arf(x,s + l) 

7+1 



3 Q 2 J 2ni V 4A 2 



(72) 



The contour integration can be computed directly after substituting Eq. (|65p or, alternatively, one can realize that 
the integral in ([72"]) is the inverse Mellin transform ([55)) at the scale /s 2 = CQ 2 /4 of the following function 



S T ^ (x. /i 2 ) J ^ rl 



" dj/ 2 
i4 



a s (p' 2 )g(x,^ 2 ) + S%' (x, nl) , 



(73) 



where the reminder (x, fly) depends only on the gluon distribution at an initial scale /Uq. It can be computed using 
the DGLAP equation ([BT3")) 



(74) 



Thus, we finally obtain 



A o- (T=2) 



4^ / r CQ2/ *d^ <2^M2) ( 2 ; 



3Q 2 



(75) 



The leading logarithmic term in Eq. ([75)) coincides with the standard DGLAP expression for ctJ p obtained assuming 
that the see quarks come from the gluon splitting in the last step of the evolution. 

The higher orders in the Laurent expansion of Ht,l(—s) m Eq. (|69p . beyond the singular term, correspond to higher 
order terms in the perturbative expansion of the twist-2 contribution. The next-to-leading order (NLO) contribution 
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originate from the constant term o!p L . The obtained expression is of the form (|72|) without (s+ 1) in the denominator. 
Thus, we immediately obtain 



A rr( T=2 ) 



NLO 3Q 



lL a s (CQ 2 /4)g(x,CQ 2 /4), (76) 



which for the transverse polarization carries one power of the large logarithm log Q 2 less than the leading term in 

(t—2) 

A a er,p . Notice that, as expected, for the longitudinal polarization the first non- vanishing twist-2 contribution 
enters at the NLO level. A similar procedure could also be applied to higher terms of the Laurent series, giving 
contributions with successively decreasing power of log Q 2 . Obviously, these higher order terms do not exhaust all the 
higher order QCD effects. They are parts of the QCD corrections to the twist-2 amplitude which come from inclusion 
of the quark transverse momentum in the quark box beyond the collinear limit 

So far we have dealt with the singular part of a(x, s) at s = —1, generated by the first term in the perturbative 
part of the Glauber-Mueller series (|58p proportional to the gluon distribution g(x,/i 2 ). The remaining terms of this 
series as well as the soft part a so ft(a;, s) are regular at s = — 1. However, they contribute to twist-2 through the pole 
of Ht(~s) at this point, giving 

A c 4^ 2) =-^5£- l) (*,* = -l). (77) 



The function on the r.h.s is a sum of two pieces: the soft part, a so {t(x,s = — 1), and the Mellin transform of the 

T (1) ( 
pert V 



regular part of the perturbative component, er pcr t(x,r) — a^^x, r), computed for ,s = —1. Thus 



) (a; ! s = -l) = J — |fr(x,r) - er^/ rt O,r)j + a soft (x, s = -1) 



(78) 



We summarize by displaying the most leading twist-2 contribution to the 7*p cross sections, obtained in the DGLAP 
improved saturation model (with C = 4): 

(r=2) 87T 2 A) 1 [ Q2 dfl' 2 , 2 , 2 

= q 752 / — 7z a s(/ i )9[x,V ) (79) 

^ =2) = ^^« s (g 2 )^,Q 2 ). (80) 

Notice the similarity concerning leading logarithms between the twist-2 contributions in the original GBW model, 
Eqs. ([50)1 and (p)Tj) . and the above formulae. 



B. Twist-4 contributions 



(2) 

The formula for twist-4 is determined by the Mellin transform &pj lt (x, s) of the second term in Eq. ([55]) and the 
Laurent expansion of Ht,s(~s) around s = —2: 

4 r 7 4) = /" ii^ S {„| ^,-)+^- a >(*,-)}{=^ + l@ ; +0(- + 2)} (81) 



2m V 4 

(i) 



where the integration contour envelopes the cut singularity with the branch point at s = — 2 in which the function 
0reg 2 ^ (x, s) is regular, see Fig. The Laurent expansion coefficient are now given by 

„( 2 ) — n ^( 2 ) — 16 /l /,(2) _ 4 4 l(2) _ (Z2 344 \ 4 /or,\ 



(2) 

The vanishing a T means that the leading logarithmic twist-4 contribution 

,(2) 



r 4 ) = oo4! r ±(or\' *) (83) 
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FIG. 6: The singularity structure in Mellin plane relevant for twist-4 transverse and longitudinal contributions together with 
the integration contours. The meaning of all symbols is as in Fig. |4j but the original contour, C 3 , is replaced by the shifted 
Cs 1 '. The fundamental strip is also shifted. 



vanishes for transverse photons. For the longitudinal polarization it can be found in a similar way as for twist-2, with 
the following result 



A J T=i ) - 8<J a L 



Q 4 V3<xo 



II 



IT, 



[ a .Qi a )g(x,n a )]' + Sf\x,nl) 



where the reminder is of non-perturbative origin and does not depend on Q 2 , 



log 



djj 



dw' 



9o(w')<?o(w-u/) 



2wi J 2m ^P gg {u>>) + ^P gg {LO - w') - 1 



(84) 



(85) 



( 21 - 

The NLO correction to twist-4 comes from the constant term, b T L , in the Laurent expansion of Ht,l{~ s ) around 
s = — 2. It is straightforward to obtain 



8ao4 2) L ( 7T 2 x J 



NLO 



Q 4 \3ctq 



[a s {CQ 2 /A)g{x,CQ 2 /A)\ 



(86) 



which in the longitudinal case has one logarithm of Q 2 less then the leading contribution (f8"4")l . As for twist-2, the 
higher terms in the Laurent expansion of Ht,l{—s) give rise to yet higher order perturbative corrections. 

Multiple scattering effects (with n > 3) and the soft contribution are important only for the longitudinal twist-4, 
crj j T_4 ' ) (x, Q 2 ). They are collected in 



16a 



(2) r 



Q 4 



d[ s ~- 2) {x, s = -2) + a soit (x, s = -2) 



where 



a r ( c s =- 2) (*> s = -2) = I ^ [a(x, r) - a« (x, r) - a< 2 > (x, r) 



(87) 



In summary, the following leading logarithmic structure is found for twist-4 (with C = 4) 

(r=4) 32tt 4 A 1 r , 2% / n 2Nl2 

4t> do y 



(89) 



(r=4) 128tt 4 A 1 f Q d\x' 2 

Tl 135ao Q 4 T 2 " 



[a s ( M /2 )3(a:,M /2 )] : 



(90) 



which should be compared to the results obtained in the original GBW saturation model, Eqs. (15"2")) and (fBU]) . Notice 
the similarity in the sign and the leading logarithmic structure. 
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C. Discussion 



The results (|79 j) . (|80|) and (|89|) .([90 |) on the leading logarithmic behaviour of the twist- 2 and twist-4 contributions 
are quite general. For the nucleon structure functions Ft and Fl they imply that the relative twist-4 correction to 
Ft is strongly suppressed w.r.t. the twist-2 contribution since the subleading twist-4 term in Ft appears only at the 
NLO. On the contrary, for Fl the leading twist term enters only at the NLO and the the twist-4 correction enters 
at the leading order. So, the relative twist-4 effects in Fl are expected to be enhanced. Note that both in the case 
of Ft and Fl, the twist-4 effects are enhanced w.r.t. the twist-2 contribution by an additional power of the gluon 
density g(x,Q ). At sufficiently small x, when the gluon density is large, this enhancement may compensate the 
twist-4 suppression w.r.t. the leading twist-2 contribution. 

For the structure function F% — Ft + Fl we expect small relative corrections from the higher twists because of the 

(2) (2) (2) (2) 

opposite sign of the terms proportional to a L and b T L . In fact, both a L and b TL are negative. Thus it follows from 
([55)) and ([50]) that the resulting LO twist-4 contribution to F2 coming from Fl is positive and both the dominant 
(though NLO) term in Ft and the NLO correction to Fl are negative. This leads to partial cancellation between the 
twist-4 LO and NLO contributions to F2 at moderate Q 2 , which can be also viewed as a partial cancellation between 
the twist-4 corrections to Fl and Ft- 

These conclusions about the importance of the higher twist corrections are expected to be quite general because 
they follow directly from the twist structure of the quark box and do not depend on the detailed form of the twist-4 
gluon distribution. In fact, for a generic twist-4 gluon density G±{x, Q 2 ) (not necessarily proportional to [g(x, Q 2 )) 2 ), 
the qualitative pattern of the computed twist-4 corrections emerges. This happens because independently of the 
detailed form of gluon density, the perturbative color dipole scattering cross section at twist-4 is proportional to 
r 4 a s (C/r 2 ) G±{x, C/r 2 ). Using a generally valid relation: f (ds/2ni) x~ s f{s)/s — f x dx' f(x'), one finds 



4 T=4) - ^a 2 s (Q 2 ) Gi (x,Q 2 ), (91) 



and 



>=4) - ' %c? s ^ 2 )G,{x^ 2 ) + ^ a 2 s (Q 2 ) Gi (x,Q 2 ). (92) 



,(2) rQ 2 a, A 
0~ J ^ 

This confirms that for twist-4 the pattern of cancellations in F2 between Fl and Ft (or between LO and NLO terms) 
is indeed universal. 



VII. TWIST DECOMPOSITION IN THE COORDINATE SPACE 



The preceding analysis was carried out in the Mellin space. This representation is useful to understand the key 
features of the twist decomposition and match the DGLAP improved saturation model with QCD. However, in the 
explicit calculations of the twist composition of the 7*2? cross sections we find it more convenient to use the coordinate 
representation. The main reason is that the multiple scattering contributions are represented as multiple convolutions 
in the Mellin space and as simple powers of r 2 a s (^ 2 )g(x, fi 2 ) in the coordinate space. Thus, we shall construct an 
explicit prescription that facilitates the twist decomposition in the coordinate space. The obtained formula (|99p is 
equivalent to its counterpart in the Mellin moment representation and it will be used to provide numerical estimates 
of the twist decomposition of the nucleon structure functions. 

The singularity structure of the product d{x, s) Ht,l{~s,Q 2 ) is similar to the structure of d(x,s) except of the 
branch points of a(x,s) which are strengthened by the poles of Ht,l(~ s iQ 2 )- I n what follows, we shall isolate the 
contributions of the singularities emerging at s = — 1, —2, . . . and link them with the twist contributions r = 2, 4, . . ., 
respectively Let us define two sets of functions, 

a^(x,r)=a(x,r) - a Q ]T [ —^n k (x,r 2 ) (93) 

k=l 

and after introducing f = rQ 

n h (k) 

45(f) = H T , L (f) E (£k ( 94 ) 
fe=i ^ ' 
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Re(s) 



FIG. 7: The singularity structure of a^°\x, s)Ht(—s,Q 2 ) and a ( - 1 \x, s)Ht{ — s, Q 2 ), and the integration contours in the 
complex s-plane. The meaning of all symbols is as in Fig. 2]- additionally, pieces of the closed integration contour are shown: 
d 1 ', rf\ etc. 



where h>p' L are the coefficients of the asymptotic expansion of Hx,l{^)i see Eq. Additionally, — a and 

H^ L — Ht,l- We see that describes the contribution to the dipole cross section of (n + 1) and more scatterings 
and Hj?' L gives the contribution to the photon wave function from the poles of twist r = 2(n + 1) and higher. The 
new functions have the following asymptotics at small and large values of r (modulo logarithms): 

and 

Now, it is easy to prove that a^- n \x,r) and iJ^(f) have Mellin transforms, a^{x, s) and H^ L {—s,Q 2 ), with the 
fundamental strip: 

- (n + 1) < Res < -n. (97) 

It is moved to the left by n units with respect to the fundamental strip given by Eq. ([24]) . 

The singularities of the Mellin transform a^ n \x, s) emerging at the branch point to the left of its fundamental strip 
are the same as the corresponding singularities of the functions: (x, s), (x, s), . . . , cr*™ -1 ) (x, s). The functions 
flj|(f) are shifted with respect to Ht,l{t) by a finite power series in 1/r 2 , so their Mellin transforms are identical 

and equal to the Mellin transform Ht,l(s, Q 2 ) for all n. The series subtraction results only in the already discussed 
shift in the position of the fundamental Mellin strip. Therefore, the difference 

A^[aH T>L ](x,Q 2 ) = / ^^-V{x,s)H£- X \- S ,C?) ~ I , a {n) (x, s) H^ ] L {-s, Q 2 ) (98) 
y c («-i) 2m J c (") 2m 

defines the contribution of the n-th singularity (i.e. the cut discontinuity with the branch point at s = —n) to the 
integral in Eq. (|2ip , see Fig. for illustration. In our analysis we identify this contribution with the r = 2n twist 

component of crj 

The Mellin integrals in Eq. (f98|) may be expressed in the coordinate space to give a formula that facilitates a direct 
determination of all twists in the coordinate representation: 

4 2 ^,Q 2 ) = J™^[o^)(x.r)H^ l \rQ) - a^(x, r) H^(rQ)} . (99) 

Clearly, the twist decomposition would be complete and 

CO 

E 4;r 2n) = 4.i (ioo) 

n=l 



(r 2 )™ 



for r 
for r 





oo 



l/f 2n for r -> 
y f 2(n+i) for r 



(95) 



(96) 
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provided that the infinite summation of A^ n > [ctHt,l\ (x, Q 2 ) is convergent 1 . 

The prescription given by Eq. (|99|) may be also applied to the original GBW dipole cross section which Mellin 
transform has a series of isolated poles at s = ~n instead of the series of cut singularities. In this case, in Eq. (|93[) a 
polynomial is subtracted and the Mellin transform of a^ n \x,r) is identical to a(x,s) given by Eq. (|^5|) . Therefore, 
formulae (|98|) and (|99|) may also be applied to extract the contribution from all singularities of &{x, s)Ht,l(—s, Q 2 ) 
to give the twist expansion in the case of the GBW dipole cross section. An explicit numerical check showed that the 
expansion obtained using prescription (|99p agrees with the analytic results in 



VIII. HEAVY QUARKS 

So far we studied the massless quark contribution. Within the fey-factorization approach it is straightforward to 
study also the case with a non-zero quark mass. In particular, the Mellin transforms of the photon wave functions 
squared with mj ^ 0, which generalize expressions (IA9|) and (|A10p . are known [12J. We shall denote them by 
Ht,l(s,Q 2 ,m 2 ). The contribution of a heavy quark to the ~f*p cross section may be obtained using the Parseval 
formula (|2ip in which the replacement Ht,l(—s,Q 2 ) — > HT,L{~s,Q 2 ,m 2 ) is made. For Res > 0, the functions 
Ht,l{s,Q 2 ,m 2 f) are regular in s. Therefore, the s-singularity structure of the integrand d(x, s)Ht.l(—s,Q 2 ,m|) in 
pip is determined by the singularity structure of s). Hence, for heavy quarks the twist-r component is determined 
by the n-fold scattering component of the dipole cross-section, 

a ( T T = 2n \x,Q 2 ,m f ) = [°° <^H T . L {r 2 ,Q 2 ,m 2 )a n {xy) (101) 
Jo r 

where a n (x,r 2 ) = <to( — l) n+1 fl n (x, r 2 )/n\. Note that for heavy quarks with m 2 ^> /it 2 ,, the integration in (| 1 1 [) does 
not lead to any infra-red divergences since the photon wave function provides an exponential cut-off proportional to 
exp(— rmf) for the r 2 integration. This was not the case for the light quarks, when m 2 <C /i 2 ,, for which formula (jlOip 
cannot be applied. 



IX. PHENOMENOLOGICAL RESULTS AND CONSEQUENCES 

In this section the obtained estimates are presented for the higher twist effects in Ft, Fl and F^, and also, separately, 
for the charm quark components of Ft , Fl and F2 . Additionally, we discuss the phenomenological consequences of our 
findings for the measurements at the LHC. We performed an explicit numerical evaluation of higher twist components 
of the proton structure functions in the DGLAP improved saturation model, defined by Eqs. l|42 p -(|43 p with the 
evolution scale given by Eq. (|54| : 

2 _ / C/r 2 for r 2 < C/fi 2 
» \ f 4 for r 2 >C/ M 2 . 

Furthermore, we compared the results to those obtained in [ll| within the GBW model [e| without the DGLAP 
evolution. 

The parameters of the DGLAP improved saturation model were fitted to describe all HERA data on F2 at x < 0.01. 
In the model we took into account three massless quark flavors and the massive charmed quark with m c = 1.3 GeV. 
The gluon density at the input scale Q 2 = 1 GeV 2 was assumed to take the form 

xg(x,Q 2 ) = A g x- x (l-x) 5 - 6 . (102) 

The parameters obtained from the best fit with \ 2 = 0.94/d.o.f are the following: 

C = 0.55, fil = 1.62, A, = 1.07, A = 0.14, cr = 22 mb. (103) 



In fact the series in not convergent; the expression was obtained assuming the validity of the asymptotic Ht l(tQ) expansion for large 
rQ down to rQ = 0. Therefore, the obtained series is asymptotic. 
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Twist ratios: tw-4/tw-2 
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FIG. 8: The ratio of twist-4 to twist-2 components of Ft, Fl and F2 at x 
lines) and in the DGLAP improved saturation model (dashed lines). 
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Structure functions 



The obtained relative twist-4 corrections (with respect to the twist-2 ones) to the structure functions Ft, Fl and 
F2 are displayed in Fig. [51 as a function of Q 2 , for x = 4 • 10~ 5 (for this value the saturation scale Q s (x) = 1 GeV 
in the GBW model with charm). The continuous curves obtained in [ll| correspond to the GBW model with charm 
quarks [T3|, and the dashed ones are obtained in the DGLAP improved saturation model (BGK) [l5[ with the 
parameters given above. The differences between the GBW model and the BGK models are visible, but rather small. 
The qualitative picture is fully consistent between the models and agrees very well with the results of the analytic 
analysis outlined in Sec. I VII Thus, the higher twist corrections are strongest in Fl, and much weaker in Ft- In F 2 
there occurs a rather fine cancellation between the twist-4 contributions to Ft and Fl, at all Q 2 , down to 1 GeV 2 . 
Although an effect of this kind was expected, it still remains somewhat surprising that this cancellation works so well. 
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FIG. 9: The ratio of the twist-2 component to the full light quark result of the model for Ft, Fl and JFa at x ■ 
GBW model with charm (continuous lines) and in the DGLAP improved saturation model (dashed lines). 
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We also show in Fig.[5]the ratio of the twist-2 component to the full dipole model result for Ft, Ft, and F 2 . The full 
result incorporates the resummed contributions of all twists. On the qualitative level, the effect of all higher twists 
shown in Fig. [9] are similar to the effect of twist-4 in Fig. [5J indicating that the higher twist corrections are driven by 
the twist-4 contribution down to Q 2 = 1 GeV 2 . It is remarkable, that the cancellation of the higher twist effects in 
F2 occurs also in the all- twist result. Clearly, the results shown in Fig. [9] indicate that at lower Q 2 , the conventional 
twist-2 calculations underestimate the value of Ft, significantly overestimate the value of Fl and slightly overestimate 
the value of F^. 

One should stress, that the theoretical conclusions about the strength of the higher twist corrections, related to the 
powers of a s , are only valid in the perturbative range, where a s is small. They are therefore, well justified for Q 2 
above, say, 5 GeV 2 . In that region, indeed, the difference between higher twists in Ft, and Ft is quite pronounced. 
At lower Q 2 , where a s is not small, the differences in powers of a s should not lead to quantitatively distinct results in 
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Ft and Fl- Indeed, at Q 2 = 1 GeV 2 , the relative twist-4 corrections to Ft and Fl are 30% and 50% correspondingly, 
that is they do not differ very much. 



B. Charm contribution 

The DIS cross-section at small- a; and a moderate Q 2 receives significant contribution from the charmed quark. The 
contribution of the bottom quark may be safely neglected due to its small charge e^ = 1/3 and its large mass. For the 
region of interest, Q 2 ~ 10 GeV 2 , the mass, m c , of the charmed quark cannot be neglected, as Q 2 ~ 4m 2 . Therefore, 
our extraction of the higher twist effects in charm structure functions F T and F£ is based on the results of Sec. IVIIII 
The results for higher twist effects in the charm structure functions F T and F£ are shown in Fig. [TQl Displayed are 
the twist- 2 and twist-4 components and the all-twist result. In contrast with the case of the light quarks, the higher 
twist effects introduce negative corrections both in F T and in F£ , and the magnitude of the ratio of twist-4 to twist-2 
contributions is similar in both cases and reaches a few percent. Consequently, the effect of higher twists in F% is 
similar. 



C. Comparison with Fl data 

Recently, new measurements were performed of the proton F L structure function in a wide kinematic range [l6| ■ The 
measurements probe F L for correlated (x, Q 2 ) pairs down to x = 5.9T0 -5 and Q 2 = 1.8 GeV 2 , see Fig. [IT] The data in 
the lowest range of x are particularly interesting, as in this region, the leading twist, fixed order DGLAP calculations 
face intrinsic problems [33. [33|. Specifically, in that region, the convergence of the subsequent F L approximations 
within perturbative expansions is rather poor, up to the next-to-next-leading order (NNLO) approximation [33j . In 
addition, at very small x and low Q 2 , the estimated Fl becomes negative, violating the fundamental condition of 
positivity [33|. This indicates that the DGLAP treatment in this region has to be improved. In what follows, we shall 
present the comparison of the new Fl data with the dipole model results, and we shall shortly compare our approach 
with another successful approach to F2 and Fl, based on the leading twist DGLAP scheme, improved by a small x 
resummation [32j. 

In Fig. [IT] we show the comparison of our results with the preliminary data on Fl from the HI collaboration. In 
the top of the plot, the values of x are indicated for each data point. Note that the experimental data points show 
strong correlation between the values of Q 2 and x. Thus, small Q 2 values are measured for smaller values of x. The 
solid curve represents the all twist result from the DGLAP improved saturation model applied in this paper, while 
the dashed line shows the twist-2 contribution within this model. The difference between the two curves comes from 
the negative higher twist terms, with a dominant contribution of twist-4. The description of the data provided by the 
model is good, both for the twist-2 approximation and the all-twist result. We stress, that all the model parameters 
are fixed by the fit to F2 data and no new parameters are introduced in the description of Fl ■ 

Clearly, the low Q 2 region of the plot, where x is also small, is highly sensitive to higher twist effects. In particular, 
for the lowest measured values of (x, Q 2 ), the twist-2 contribution is already about two times larger then the exact 
result. Unfortunately, the current experimental errors are sizable and no evidence for higher twist effects can be drawn 
from the measurements, yet. We stress, however, that Fl at small x and Q 2 should be an excellent observable to find 
such effects, provided that the experimental errors may be further reduced. 

The defects of the fixed order DGLAP description of Fl at small x and Q 2 were shown to be partially cured 
by including into the DGLAP framework a resummation of small x corrections, enhanced by powers of logs. The 
resummation, proposed by Thorne and White (TW) [321] . absorbs the NLL BFKL effects at the leading twist into the 
NLO DGLAP evolution. In the currently relevant kinematic range, the description of Fl based on the TW approach 
is remarkably similar [331 to the one obtained within a saturation model with the impact parameter dependence (the, 
so called, b-Sat model) [3J|. In addition, the b-Sat model results for Fl agree well with the results of this paper. 
The TW scheme provides a good description of the existing Fl data. The Fl at small x and Q 2 following from the 
saturation models is significantly lower than the corresponding TW result, but the differences are not pronounced 
enough to permit a discrimination between the approaches with the present data. Let us, however, stress, that the 
asymptotic x — ► (or Q 2 — > 0) behaviour of the structure functions should be different in approaches consistent with 
unitarity constraints, (as e.g. the saturation models), and the leading twist approach. In the former case, Fl should 
vanish in the limit, while in the latter case it should remain non-zero. Thus, one expects, that the leading twist 
approach should be insufficient at a very low x and fixed Q 2 , and that the inclusion of higher twist effects should be 
necessary in that limit. 
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FIG. 10: Charm contribution to Ft and Fl at x — 4 ■ 10 -5 in the DGLAP improved saturation model. Shown are the 
contributions of the twist-2 component (dashed lines), the twist-4 component (dotted lines) and the full result (continuous 
lines). 



D. Discussion and implications for the LHC 



The analysis performed in this paper shows that the importance of the higher twist corrections may essentially 
depend on the process and the probe. In particular, higher twist effect in the structure function F2 are strongly 
suppressed due to rather fine cancellations between Ft and Ft,. Such cancellations are not expected to occur in a 
generic case. For instance, the higher twist effect in Fl are enhanced. Thus far, parton density functions (pdfs) in 
DIS were fitted mostly to the F2 data. Due to small higher twist effects in F%, one expects that the suppressed higher 
twist contributions should not affect the quality of the determination of pdfs. This is, certainly, a good news. The 
estimated correction due to higher twist effects in Fl at small x and a moderate Q 2 is, however, much larger, and 
this correction should be taken into account when including the Fl data into fits of pdfs. 

The example of F2 and Fl in the DIS shows that the multiple scattering (higher twist) effects are probed in various 
ways, depending on the observable. Similar differences in the magnitude of higher twist effects in various observables 
may occur in the hadronic collisions, e.g. in pp collisions at the LHC. In particular, cancellations present in F2 is not 
expected for the bulk of LHC observables probing the gluon distributions at small- a;. Thus, in general, the relative 
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Comparison with H1 data 
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Q 2 /GeV 2 

FIG. 11: Comparison of the DGLAP improved saturation model with the preliminary HI on Fl [Hj]. The solid line is the 
all- twist result while the dashed line shows the twist-2 contribution. The value of x is indicated for each data point. 

effects of higher twists at the LHC should be larger than they are in F 2 ■ As an example, let us give the case of the 
forward Drell-Yan process, that can be effectively described using the dipole formulation [35]. At LHC-b, the Drell- 
Yan process may be probed at moderate Q 2 ~ 10 GeV 2 and x ~ 10 -6 , what should provide a gold-plated probe of 
the gluon density at small x. However, the higher twists effects may be quite strong there. In particular, let us invoke 
an example of the Lam-Tung relation [3(| that holds for angular distribution of Drell-Yan lepton pairs. According to 
this relation, the twist-2 contribution to one of the angular components of the dilepton distribution vanishes in the 
leading logarithmic approximation. Therefore, higher twist effects in this component should be enhanced, in analogy 
to the case Fl [36l. [37j. Thus, given the low values of x and Q 2 , that can be reached in the measurements in LHC-b, 
the violations of Lam- Tung relation should provide a sensitive probe of the higher twist gluonic operators at small x. 
On the other hand, the higher twist effect may be also large in the total cross-section of the forward Drell-Yan process. 
In that case, a determination of gluon density at small x, based on the leading twist contribution alone would be 
inaccurate, and the higher twist contributions should be taken into account. Besides that, the higher twist effects 
may be larger in processes with gluons, like e.g. the forward gluonic jet productions, where the multiple scattering of 
the gluon is enhanced by its color charge, as compared to the quark rescattering in the DIS case. In such processes, 
we do not expect that any cancellations of rescattering effects should occur, of the type found in F 2 . 



X. CONCLUSIONS 

In this paper the leading higher twist contributions to proton structure functions, F2, Ft and Fl, at small Bjorken x 
and moderate Q 2 were analyzed. The problem was analyzed theoretically confronting two different approaches. In 
the fist approach, we focused on a subset of QCD diagrams describing contributions of quasipartonic gluon operators, 
that should dominate the higher twist effects in the deeply inelastic scattering at small x. We demonstrated, that this 
subclass of the diagrams, at the leading logarithmic approximations and in the large N c leads to a picture consistent 
with the DGLAP improved saturation model. In contrast, we considered also the problem of higher twists in the 
Balitsky-Kovchegov framework, in which, the BFKL Pomeron fan diagram are resummcd. In this approach, higher 
twist contributions coming from the fan diagrams vanish in the leading log Q 2 approximation. 

The pattern of the most important twist-2 and twist-4 contributions to F%, Ft and Fl is determined by the properties 
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of the quark loop through which the virtual photon interacts with the gluonic field of the proton. Therefore, it is 
universal and its key features should not depend on the model details. Those features are: (i) the twist-4 correction 
to Ft enters only at the NLO, and so, the twist-4 correction to Ft is suppressed; (ii) the twist- 2 contribution to Fl 
enters at NLO, and the LO twist-4 term in Fl is relatively enhanced and more important; (iii) the relative sign of 
twist-4 corrections to Ft and Fl is opposite, and the higher twist effects partially cancel in F2 = Ft + Fl- These 
general conclusions were then confirmed by a quantitative phcnomenological analysis. 

We performed a numerical twist analysis of the DIS cross-sections at small- x within the DGLAP improved saturation 
model. In order to carry out a quantitative estimate of the higher twist effects in the structure functions, we proposed 
a method allowing for a direct, numerical twist decomposition of the saturation model cross sections. The method 
was then applied to the DGLAP improved saturation model, fitted to the HERA F2 data. Contributions of twist-2, 
twist-4 and all twists to F2, Ft and Fl were then extracted. We found a strikingly good cancellation of the higher 
twist effects in F2, for which, at x = 3 • ICP 4 , the relative correction from higher twists is found to be at a few 
percent level down to Q 2 — 1 GeV 2 . The higher twist corrections to Ft were found to be moderate, below 10% for 
Q 2 > 3 GeV 2 at x = 3 • 1CP 4 . On the other hand, the twist-4 correction in Fl was found to be large, about 50%, at 
Q 2 = 1 GeV 2 , and still sizeable, about 20%, at Q 2 = 10 GeV 2 . Therefore, whereas the leading twist analyses of F 2 
are fully justified, one should include the higher twist effects in analyses of the Fl data at small x and moderate Q 2 . 
We also found that the saturation model description of the recent Fl measurements at small x and low Q 2 is good. 
Unfortunately, the data are not precise enough to prove that the inclusion of higher twist corrections improves the 
description of the data. 

Finally, some implications were discussed of the results for analyzes of the LHC data. In particular, we stressed 
a strong process-dependence of the higher twist contributions, exemplified before by the striking differences between 
F2 and Fl- It follows from our analysis, that F2 is protected by cancellations from the higher twist effects, and 
such cancellations are not expected to be generic. Therefore, the higher twist effects in some LHC observables may 
be much stronger that they are in F2. Thus, it is crucial to estimate higher twist effects when attempting a precise 
determination of parton densities in LHC measurements at small x and moderate Q 2 , like e.g. in the forward Drell-Yan 
process at low Q 2 , or in the forward jet production. 
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APPENDIX A: MELLIN TRANSFORMS OF H t ,l 

Let us compute the Mellin transform of Ht(tQ) given by Eq. (|23|) 

/>OC 

H T (s,Q 2 ) = dr 2 (r 2 )*- 1 H T (rQ) . (Al) 



Jo 

Substituting Eq. ([23]) we obtain 

v 4 

where 

1 



Ht{s,Q 2 ) =(^1 II r(,) (A2) 



H T (s)=4- s A dz [z 2 + (I - z) 2 } z{l - z) dr 2 (f 2 ) s K\(sJ z(\ - z) f) . (A3) 
Jo Jo 

with r = rQ. Changing the variable, y 2 = z(l — z)f 2 , we find 

r 1 z 2 + n - z) 2 r°° 

H T ( S )=4- S A dz—L} J- dy 2 {y 2 YK 2 {y). (A4) 

Jo z s {l-z) s J 

The integral over z equals 

Jo z*(i- z y J ^ ; r(4-2 S ) 



(A5) 
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where we used the definition of the Euler beta function. The gamma function in the denominator can be written as 

o2(2-s)-l 



r(2(2- s )) = 



• r(2 - s)r(2 - s + 1/2) 



and from this we have 



_ r(i- a )r(3- a ) 

z A 1 -" r(2 - s)I\5/2 - s) 



The integral over y in Eq. (|A4p equals 



Thus, we finally find 



- r(2 + a) r(i + a) r(«) r(i - s ) r(3 - s) 

H T (s) - 



8 r(3/2 + s)T(2- s)T(5/2- s) 
A similar calculation allows to compute the Mellin transform 

a q tt (r(i + s )) 3 r(2- s ) 



4 r(3/2 + s)T(5/2- s) 



(A6) 



(A7) 



(A8) 



(A9) 



(A10) 



APPENDIX B: EVOLUTION OF THE GLUON DENSITY 



The gluon density used in our analysis, g(x,p 2 ) = xG(x,fi 2 ), obeys the following leading logarithmic DGLAP 
evolution equation: 



dp? 2?r j x \ z 

where the contribution from quarks is neglected and the gluon splitting function P gg takes the form, 



PM = 6 



+ z(l -z) + - + -6(1 - z) 



and the one loop strong coupling is given by 



a s (M 2 ) 



(l-z)+ 12 



b 



log(M 2 /A 2 ) 

with bo = 127r/(33 — 2ri/). This equation may be solved in the double Mellin moment space, 



where g(uj,j) obeys the following equation 



dio 



2^ 5(W ' 7) U 



2 \ 7 



^;{73(^,7)} = ^ S9 (w)j(u,7) 



and the splitting kernel in the Mellin representation, P gg (uj), is given by 



P gg (o>)= / dzz"P gg (z) = 6 



1 1 



lo u) + 1 u> + 2 uj + 3 



j E - f/>(u + 2) 



33 - 2n f 
6 



The general solution to Eq. (|B5[) reads 



(Bl) 



(B2) 



(B3) 



(B4) 



(B5) 



(B6) 



(B7) 
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where go(uj) is an arbitrary function which may be fixed using an initial condition. Thus, the solution expressed in 
terms of the original variables (x, jjl ) is given by 



djj 



9( x ^ ) = / -K- x "9oH / 7—1 



2m 



2m 



,2\ 1 



A 2 



(B8) 



The contour integral over 7 may be performed for all /i 2 > A 2 after the integration contour is deformed to envelope 
the cut along the negative real half-axis. We find 



du> 
2~Ki 



x g{w,ii ) 



(B9) 



where 



<7oM 



log 



(BIO) 



The initial condition for the DGLAP equation at some scale (1q ^> A 2 is given by its Mellin transform g(u>, /i 2 ,). Thus, 
writing (jBlOj) for it 2 = /z 2 ,, we obtain 



log 



(Bll) 



which in turn, after the substitution to (|B9|) . leads to the well known form 



log( M 2 /A 2 ) 



log(Ai 2 /A 2 ) 



(B12) 



We also need the double Mellin representation of the product a s {^ 2 )g(x^i 2 ). In the mixed (uj,fi 2 ) representation 
the DGLAP equation (jBl~j) reads 



2 dg(uj,H 2 ) p gg(u) , 2 x -/ 2n 

M — tt^ = — ^ <W 

ctyr 27r 

Taking the Mellin moment (|6"2"|) of both sides (with s = 7), we obtain 



(B13) 



7<?0>7) 



Pgg{") 
2tt 



a s g{uJ,l) ■ 



(B14) 



Thus, after inserting relation (|B7j) we find 



«s5 



. . 2m]Q(uj) i> p , . 



which after coming back to the x variable reads 



(B15) 



u s g{xn) 



2™ p m M 



(Bi6) 



Both functions have logarithmic cut singularity in the complex 7-plane with the branch point at 7 = 0. 
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